home *** CD-ROM | disk | FTP | other *** search
/ EnigmA Amiga Run 1995 October / EnigmA AMIGA RUN 01 (1995)(G.R. Edizioni)(IT)[!][issue 1995-10][Aminet 7].iso / Aminet / dev / gcc / ixemul_src.lha / ixemul-41.0 / stdlib / atof.c < prev    next >
C/C++ Source or Header  |  1995-05-17  |  3KB  |  122 lines
  1. /*-
  2.  * Copyright (c) 1990 The Regents of the University of California.
  3.  * All rights reserved.
  4.  *
  5.  * This code is derived from software contributed to Berkeley by
  6.  * the Systems Programming Group of the University of Utah Computer
  7.  * Science Department.
  8.  *
  9.  * Redistribution and use in source and binary forms are permitted
  10.  * provided that: (1) source distributions retain this entire copyright
  11.  * notice and comment, and (2) distributions including binaries display
  12.  * the following acknowledgement:  ``This product includes software
  13.  * developed by the University of California, Berkeley and its contributors''
  14.  * in the documentation or other materials provided with the distribution
  15.  * and in all advertising materials mentioning features or use of this
  16.  * software. Neither the name of the University nor the names of its
  17.  * contributors may be used to endorse or promote products derived
  18.  * from this software without specific prior written permission.
  19.  * THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY EXPRESS OR
  20.  * IMPLIED WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED
  21.  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
  22.  */
  23.  
  24. #if defined(LIBC_SCCS) && !defined(lint)
  25. static char sccsid[] = "@(#)atof.c    5.1 (Berkeley) 5/12/90";
  26. #endif /* LIBC_SCCS and not lint */
  27.  
  28. #define KERNEL
  29. #include "ixemul.h"
  30.  
  31. #include <ctype.h>
  32.  
  33. #ifdef amigados
  34. #define IEEE
  35. /* seems like the amiga stores IEEE in the same byte order as the 300, 
  36.  * lucky we:-)) */
  37. #define hp300
  38. #endif
  39.  
  40. static const double _twoemax =
  41. #ifdef IEEE
  42.     9007199254740992.;    /*2^53*/
  43. #else
  44.     72057594037927936.;    /*2^56*/
  45. #endif
  46.  
  47. #ifdef hp300
  48. /* attempt to be as exact as possible */
  49. static const struct {
  50.     long d_high;
  51.     long d_low;
  52. } _exp5[] = {
  53.     { 0x40140000, 0x00000000 },    /* 5 */
  54.     { 0x40390000, 0x00000000 },    /* 25 */
  55.     { 0x40838800, 0x00000000 },    /* 625 */
  56.     { 0x4117d784, 0x00000000 },    /* 390625 */
  57.     { 0x4241c379, 0x37e08000 },    /* 152587890625 */
  58.     { 0x4493b8b5, 0xb5056e17 },    /* 2.3283064365387e+022 */
  59.     { 0x49384f03, 0xe93ff9f6 },    /* 5.42101086242753e+044 */
  60.     { 0x52827748, 0xf9301d33 },    /* 2.93873587705572e+089 */
  61.     { 0x65154fdd, 0x7f73bf3f }    /* 8.63616855509445e+178 */
  62. };
  63. #else
  64. static const double    _exp5[]    = {
  65.     5.,
  66.     25.,
  67.     625.,
  68.     390625.,
  69.     152587890625.,
  70.     23283064365386962890625.,
  71. #ifdef IEEE
  72.     5.4210108624275231e+044,
  73.     2.9387358770557196e+089,
  74.     8.6361685550944492e+178,
  75. #endif
  76. };
  77. #endif
  78.  
  79. double
  80. atof (const char *p)
  81. {
  82.     extern double ldexp();
  83.     register c, exp = 0, eexp = 0;
  84.     double fl = 0, flexp = 1.0;
  85.     int bexp, neg = 1, negexp = 1;
  86.  
  87.     while((c = *p++) == ' ');
  88.     if (c == '-') neg = -1;    else if (c == '+'); else --p;
  89.  
  90.     while ((c = *p++), isdigit(c))
  91.         if (fl < _twoemax) fl = 10*fl + (c-'0'); else exp++;
  92.     if (c == '.')
  93.     while ((c = *p++), isdigit(c))
  94.         if (fl < _twoemax)
  95.         {
  96.             fl = 10*fl + (c-'0');
  97.             exp--;
  98.         }
  99.     if ((c == 'E') || (c == 'e'))
  100.     {
  101.         if ((c= *p++) == '+'); else if (c=='-') negexp = -1; else --p;
  102.         while ((c = *p++), isdigit(c)) eexp = 10*eexp + (c-'0');
  103.         if (negexp < 0) eexp = -eexp; exp += eexp;
  104.     }
  105.     bexp = exp;
  106.     if (exp < 0) exp = -exp;
  107.  
  108.     for (c = 0; c < sizeof(_exp5)/sizeof(_exp5[0]); c++)
  109.     {
  110. #ifdef hp300
  111.         if (exp & 01) flexp *= *(double *)&_exp5[c];
  112. #else
  113.         if (exp & 01) flexp *= _exp5[c];
  114. #endif
  115.         exp >>= 1; if (exp == 0) break;
  116.     }
  117.  
  118.     if (bexp < 0) fl /= flexp; else fl *= flexp;
  119.     fl = ldexp(fl, bexp);
  120.     if (neg < 0) return(-fl); else return(fl);
  121. }
  122.